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j—j ' Abstract. Studies of the mass function (MF) of open clusters of different ages allow 

us to probe the efficiency with which brown dwarfs (BDs) are evaporated from clusters 
to populate the field. Surveys in old clusters (age > lOOMyr) do not suffer so severely 
from several problems encountered in young clusters, such as intra-cluster extinction 
and large uncertainties in BD models. Here we present the results of a deep photo- 
metric survey to study the MF of the old open cluster Praesepe (age 590^20 Myr and 
distance 190*j9pc), down to a 5<x detection limit at i ~25.6 mag (~40Mj up ). We iden- 
i tify 62 cluster member candidates, of which 40 are substellar, from comparison with 

predictions from a dusty atmosphere model. The MF rises from the substellar boundary 
until ~60 Mj U p and then declines. This is quite different from the form inferred for other 
f^*) ■ open clusters older than 50 Myr, but seems to be similar to those found in very young 

open cluster, whose MFs peak at ~10Mj up . Either Praesepe really does have a different 
MF from other clusters or they had similar initial MFs but have differed in their dynam- 
ical evolution. We further have identified six foreground T dwarf candidates towards 
Praesepe, which require follow-up spectroscopy to confirm their nature. 

i 

1. Introduction 



The mass functions (MFs) of stellar and substellar populations have been determined 
from optical and near-infrared surveys for several open clusters at different ages, such 
as the Orion Nebula Cluster, cr Orionis, p Ophiuchi, Taurus, IC 348, IC 2391, M35, 
the Pleiades, and the Hyades. Studies of relatively old open clusters (age > 100 Myr) 
are important for the following two reasons in particular: first, they allow us to study 
the intrinsic evolution of basic properties of BDs, e.g. luminosity and effective temper- 
ature, and to compare the evolution with structural and atmospheric models; second, 
we may investigate how the BD and low-mass star populations as a whole evolve, e.g. 
the efficiency with which BDs and low-mass stars evaporate from clusters. Such an 
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investigation has been car ried out for the Hyade s ( Bouvier et al.ll2008l and references 
therei n) and for Praesepe dBoudreault et alJl2O10L and references therein). 

iBoudreault et al.1 (120 id . hereafter B2010) observed a significant difference be- 
tween the MFs of Praesep e and Hyades: Wh ile the Hyades MF is observed to have 
a maximum at -0.6 M Q dBouvier et all 120081) . the MF of Praesepe continues to rise 
from 0.8 M Q down to 0. 1 M©. This is surprising, as both clusters share similar physical 
properties (ages, mass, metallicity, and tidal radii). Disagreement between the Praesepe 
and Hyades MFs could arise from varia tions in the cluste rs' initial MFs, or from dif- 
ferences in their dynamical evolution (IBastian et al.ll2010h . Although different binary 
fractions could cause the observed (system) MFs to differ, there is no clear evidence for 
varying binary fractions from measurements published in the literature (B2010). 



2. Observations and analysis 

The Large Binocular Cameras (LBCs) are two wide-field, high-throughput imaging 
cameras, namely Blue (LBCB) and Red (LBCR), located at the prime focus stations of 
the Large Binocular Telescope (LBT). Each LBC has a wide field of view (23'x23'), 
with four CCD detectors of 2048x4608 pixels each, providing images with a sampling 
of 0.23"/pixel. The optical design and detectors of the two cameras are optimized for 
different wavelength ranges: one for ultraviolet-blue wavelengths (including the Bessel 
U, B, V and Sloan g and r bands), and one for the red-infrared bands (including the 
Sloan i, z and Fan Y bands). In the full binocular configuration, both cameras are 
available simultaneously, and both point in the same direction of the sky, thus doubling 
the net efficiency of the LBT. The survey was carried out with the r filter using LBC- 
blue and the iz,Y filters using LBC-red, covering the central 0.59 deg 2 area of Praesepe. 

The standard data reduction steps for the LBT data were performed using the IDL 
astronomy package and IRAF. An astrometric solution was achieved using the Sloan 
Digital Sky Survey (SDSS) catalogue as a reference. The root mean square accuracy of 
our astrometric solution is 0.10-0.15 arcsec. To correct for Earth atmospheric absorp- 
tion on the photometry, we calibrated the infrared data using the r, i and z band values 
of SDSS objects which were observed in the science fields. In order to calibrate our 
Y band photometry, we used our LBT i and z photometry and the Y band photometry 
from the United Kingdom Infrared Telescope Infrared Deep Sky Survey (UKIDSS). 



Candidate Selection Procedure and mass determination 



The c andidate selection procedure and the mass determination introduced by lBoudreault & Bailer- Jonesl 
(2009) and B2 010 were ad opted in the present work. We use t he evolutionary tr acks 
from dChabrier et alj |2~000) and the atmosphere models from ( Allard et alj|200lb as- 



suming a dusty atmosphere (the AMES-Dusty model), to compute an isochrone for 

^° Myr, a distance of 190+f ° ] 



Praesepe using an age of 590^™ Myr, a distance of 190^*2 pc, a solar metallicity and 



we neglecting the reddening. 

3.1. Candidate Selection using colour-magnitude and colour-colour diagrams 

Candidates were first selected from our CMD by keeping all objects which are no more 
than 0.28 mag redder or bluer than the isochrones in all CMDs. This number accom- 
modates errors in the magnitudes and uncertainties in the model isochrones. We also 
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Figure 1 . Colour-magnitude diagram (Left ) and colour-colour diagram (Right) 
used for the first and second selection procedures. Solid lines are the isochrones 
computed from an evolutionary model with a dusty atmosphere (AMES-Dusty). The 
dashed lines delimit our selection band. The numbers indicate the masses (in Mj up ) 
on the model sequence. In the right panel, the theoretical colours of six galaxies and 
of red giants are shown as thin lines and as thick lines, respectively. The six galaxies 
are two starbursts, one Sab, one Sbc, and two ellipticals of 5.5 and 15Gyr, with 
redshifts from z-Q to z=2 in steps of 0.25 (evolution not considered). We assume 
that all red giants have a mass of 5M Q , 0.5 < log g < 2.5 and 2000 K < T e g < 
6000 K. 



include the errors from age and distance of Praesepe. We additionally include objects 
brighter than the isochrones by 0.753 mag in order to include unresolved binaries. In 
Figure Q] (left) we show the CMD where candidates were selected based on z vs. i-z. 

The second stage of candidate selection involves retaining just those objects which 
lie within 0.28 mag of the isochrone in the colour-colour diagram. This value reflects 
the photometric errors and uncertainties in the model isochrones. In addition, we also 
include the uncertainty in the age estimation of Praesepe. The colour-colour diagram 
with the selection limits is shown in Figure[T](right ), with the th eoretical colours for red 
giants using the atmosphere models of Hauschildt et al. ( 1999) and theoretical colours 
of six galaxies with redshift from to 2 from K. Meisenheimer et al. ( in prep. ) over- 
plotted. Neither the red giants nor the galaxies are expected to be a significant source of 
contamination; most of the low redshift galaxies can be easily rejected through visual 
inspection. As these are the dominant potential contaminators, we conclude that there 
is no significant contamination of non-Praesepe members in our sample (Praesepe is at 
a Galactic latitude of b = +32.5°). 

3.2. Observed magnitude vs. predicted magnitude 

Our determination of T e fi is based on the spectral energy distribution of each object 
and is independent of the assumed distance. The membership status of an object can 
therefore be assessed by comparing its observed magnitude in a band with its magnitude 
predicted from its r e ff and the Praesepe's isochrone (which assumes a distance and an 
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age). This selection step is only a verification of the consistency between the physical 
parameters obtained of the photometric cluster candidates with the physical properties 
assume for the cluster itself when computing the isochrones. In order to avoid removing 
unresolved binaries that are real members of the cluster, we keep all objects with a 
computed magnitude of up to 0.753 mag brighter than the observed magnitude. This 
selection procedure is illustrated in Fig. El 
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Figure 2. Difference between the observed J magnitude and the model J magni- 
tude computed from the derived mass and T e g, as a function of T eJ f. The vertical line 
marks the location of L0 dwarfs. The dotted line (at -0.753 mag) represents the off- 
set due to the presence of unresolved binaries, the dashed-dotted lines represent the 
error on the magnitude determination, and the long-dashed lines for the uncertainties 
on the age and distance of Praesepe. The horizontal solid line just traces zero. 



4. Results 



4.1. Selected photometric candidates 



62 photometric candidates survive the selection procedures (based on isochrones as- 
suming dusty atmospheres). This corresponds to ~ 105 objects per square degree. Our 
survey saturation occurs at ~ 18 mag in z band, corresponding to ~100 Mj ]ffi . There- 



fore, most of the low mass candidates discovered in previous surveys (e.g. Pinfieldetal 



i997LlHamblvetal.lll995h saturate in our LBT images . Only a few faint BDs classified 



by iPinfield et al. I (ll997h . lGonzalez-Garcfa et all (l2006h and by B2010 are rediscovered 



by the current survey. 

There are also some targets which are previously identified as cluster members but 
rejected by our selection procedures or visual inspection. For example, ten candidates 
identified by B2010 are detected (not saturated) in our LBT survey, but six of them are 
rejected by the z vs. i-z CMD, because they are bluer than the isochrones area. Another 
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one is obviously not a point-like source in the LBT image, and another is rejected 
because its observed J magnitude is not consistent with its model-predicted magnitude. 
Only the remaining two targets are confirmed to be cluster dwarf stars. As the current 
work employed more photometric bands than B2010 did, it is not surprising that we 
achieve a more conservative selection. 

Most of our candidates are in the substellar regime, and other than these ten, no 
other accurate photometric observations are available from past epochs. This precludes 
using proper motions as a mean of cluster membership assessment at this time. 

4.2. Contamination by non-members 

As mentioned before, the two main sources of contamination are the background red 
giants and unresolved galaxies. Red giants occupy the high mass end of this study, as 
seen in the i - J vs. z - K s diagram. Although some types of galaxies share similar 
colours with Praesepe cluster members more massive than 60Mj up , such low-redshift 
galaxies are in general extended sources and therefore easily rejected by our visual 
inspection. Among the 74 candidates that passed our selection procedures, we identify 
four as galaxies through their LBT images. Other possible sources are field L dwarfs 
and high redshift quasars (for instance at z ~6: ICaballero et al1 l2008). However, as such 
quasars have spectral energy distributions similar to mid-T dwarfs whereas our faintest 
candidates have colours of early L d warfs, and given t hat they are rare (0.25 quasars at 



5.5 < z < 6.5 in a 0.59 deg 2 survey, iStern et al.H2007l) . the MF should not be affected 
by quasar contamination. 



The contamination by L dwarfs is also unimportant. ICaballero et al.1 (2008) have 



collected possible field dwarf contaminants covering spectral type from M3 to T8 from 
the literature. From their Table 3, the spatial density for L dwarfs in the solar neigh- 
bourhood is ~ 7 x 10~ 3 pc~ 3 . Given that Praesepe has a Galactic latitude of +32.5 deg 
and distance of 190 pc, the nearby spatial density of L dwarfs Praesepe should be 
~ 6 x 10~ 3 pc~ 3 , assuming an exponential decrease for stellar density perpendicular 
to the Galactic disk and a scale height of 500 pc. If we define a volume corresponding 
the area of our survey, and use the distance uncertainties to the cluster as its depth, the 
total volume is ~ 80 pc 3 . Therefore, we estimate that we have ~0.5 L dwarf contami- 
nants near the cluster, which amounts to a negligible contamination of merely 0.7%. A 
similar calculation shows that we would have ~4 M dwarf contaminants, about 6%. 

We conclude that various contaminants are not important for this study and the 
MF we derive for Praesepe should be accurate. 



5. Mass function of very low mass and substellar population of Praesepe 

The mass function, ^(logioM), we present here is the total number of objects per square 
degree in each logarithmic mass interval logioM to logioM + 0. 1, Since we do not make 
any corrections for binaries, we compute here a system MF. 

Our optical photometry reaches lower masses than the NIR photometry that we 
used. To compute the MF of Praesepe to the lowest mass bin, we first computed a MF 
using only the optical iz photometry. This MF is presented on Fig. [3] as filled dots. 
We computed a second MF from the list of candidates that pass the three selections 
criteria which are also detected in the survey of B2010 in the NIR J and K s bands 
(presented on Fig.|3]as filled triangles). For each mass bin, we computed the number of 
object removed because of adding the J and K s filters to our selection process and mass 
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determination (plotted as a function of mass in Fig. [3j top panel). We fitted a linear 
function to estimate the number of objects that would be removed if we had additional 
/ and K s photometry down to 40^15 Mj up . The corresponding extension of the MF with 
izJK s photometry is given as a large triangle on Fig. [3] 
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Figure 3. Bottom. MFs based (a) on our survey LBT i and z photometry (dots), 
and (b) also combined with the J and K s photometry from B2010 (triangles). Error 
bars are Poissonian arising from the number of objects in each bin, except for the 
last bin, for which the error bar is mostly from the linear fit in top panel. The vertical 
dotted line is the saturation mass. The vertical long and short dashed lines are the 
masses at the 5cr detection limits of our optical LBT data and of the B2010 NIR data. 
Top. Difference of the number of objects in each mass bin, between the MF computed 
using the LBT iz data and the MF computed using the combination of the iz data from 
the B2010 NIR JK S data. The dotted line is a linear fit to the discrepancies. 



Our derived MF of Praesepe (presented in Fig. 3) shows a rise from 105 Mj up 
to 60Mj up and then a turn-over at ~60Mj. We note that in the second mass bin at 
~ 48 Mj up , the MF is very low. This is probably because at this mass we reach the 5cr 
detection limit of the J-band of our £l2k photometry, so we suspect this is not a real 
feature. However, the turn-over at ~60 Mj up occurs well above the 5cr of either iz bands 
or JK S bands (e.g., at ~60Mj up , i ~ 22, completenss~ 100%) and therefore should 
not be caused by incompleteness. This is the first time one observes a clear rise in the 
substellar MF in old open cluster. 

We collected some re sults for other clusters for compa risons, as plotted in Fig. |4l 

This i ncludes: IC 2391 by dBoudreault & Bailer- Jonesi r2009h. ON C by (iHillenbrand & Carpenter 
20001) . o~ Orionis by (Bihain_etaL||2009j) and the Hyades from (iBouvier et al.ll2008h . 
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Figure 4. Mass functions of various open cluster. From top to bottom: pOph 
(~lMyr); S 106 (~ 1 Myr, middle panel for example); Trapezium (-0.8 Myr); 
IC2391 (-50 Myr); ONC (~5 Myr); <x Ori (~3 Myr); Hyades (-625 Myr); Praesepe 
(-590 Myr, from this work). We also show the Galactic field star MF from 
(2003) as a dashed line. All the MFs are normalized to the lognormal fit of 
(12003b at the substellar limit (-72 M Jup ). 
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The MF of Praesepe is quite different from either IC 2391 ( age of ~ 50 Myr) or the 
Hyades (-625 Myr). Either the 'dynamical evaporation' does not have the same effect 
on those three clusters, or they have a different initial mass function. Another alterna- 
tive possibility is that Praesepe has a different binary fraction. Further studies would be 
necessary to clarify these points. 

The continuing rise of the MF into the substellar regime which we observe has 
also been observed in young cl usters (as shown in Fig. @}, specifically in <r Orionis 
(Biha in et alj2009h . Trapezium (|Muench et al.ll2002l turn-over at -10-20 M Jup ), p Oph 
dMarsh et al.ll2010l. MF rising to -10 Mj) and in the very low luminosity yo ung cluster 
in S1 06, where the MF increases or at least remains flat down to ~10Mj up ( Oas a et al.1 



2006). If we assumed a universal IMF, then it seems that the substellar MF of Praesepe 



has not evolve significantly since it the cluster formed. 



6. Conclusions 

We have carried out the deepest survey to date of the open cluster Praesepe, covering 
the central 0.59 deg 2 in the rizY bands. The survey probed a mass range from -100 to 
40Mj up at 5cr detection limit. Our /z-bands data, combined with the QZk NIR (7 & 7^ s ) 
band observations from B2010, are compared with theoretical loci of cluster members 
based on a dusty atmosphere (the AMES-Dusty model), in order to clarify cluster mem- 
ber candidates. Our final sample comprises 62 photometric candidates. We estimate the 
contamination by field L dwarfs to be less than 1%, and that by galaxies and red giants 
also to be negligible. About two thirds of our candidates have theoretical masses below 
the Hydrogen-burning limit at 0.072 M Q , and are therefore BD candidates. 
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The mass function we have inferred for Praesepe is consistent with that inferred 
by B2010 at a mass just below the substellar boundary, but deviates by ~0.4dex in the 
next lowest mass bin, which may indicate either a significant number of objects missing 
in the Boudreault et al. 2010 survey, or a higher concentration of substellar objects in 
the centre of Praesepe (as the Boudreault et al. survey is at a larger cluster radius). The 
latter possibility suggests that the dynamical evolution of very low mass stars is not 
efficient in this cluster, as proposed by B2010, for explaining the discrepancy between 
the Praesepe MF and Hyades MR 

The steady rise of the Praesepe MF until ~60Mj up was unexpected. Such a signif- 
icant peak has never been observed in any other cluster older than 50 Myr, but has been 
observed in several very young open clusters such as cr Orionis or clusters in star form- 
ing regions (e.g. Trapezium). This suggests that the dynamical interactions in Praesepe 
have very little effect on MFs, if we assume there is a universal initial MF. 

The results reported here will be presented with further details in an future publi- 
cation submitted to A&A (Wang et al. 2010, submitted). 

Acknowledgments. Some of the observations on which this work is based were 
obtained during LBT programme "LBT-F08-02". 
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